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Wc search for B meson decays into two-body combinations of q, r\ , uj, and (f) mesons from 89 
million BB pairs collected with the BABAR detector at the PEP-II asymmetric-energy e'^e~ collider 
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at SLAC. We find the branching fraction B{B° r]uj) = {4:.0ti^ ± 0.4) x 10"'^ with a significance 
of 4.3 a. For the other decay modes we set the following 90% confidence level upper limits on the 
branching fractions, in units of 10"^ B{B° -> rirj) < 2.8, B{B° -» rirj') < 4.6, B{B° -> r]'ri') < 10, 
B{B° ri'to) < 2.8, B{B° ^ < 1.0, B{B° -q' 4>) < 4.5, and B{B° ^ <t)4>) < 1.5. 



PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 

We report the results of searches for B'^ meson decays 
to two charmless pseudoscalar mesons ry?/, 7777', ri'rj' , to 
the pseudoscalar- vector combinations 7]uj^ 77' w, tjc^, 77' 0, 
and to the vector meson pair (jxp. These together with 
LJLO and LO(j) constitute all combinations involving isospin 
singlet members of the ground state pseudoscalar and 
vector-meson nonets. These decay modes have not been 
observed previously; the published experimental upper 
limits on their branching fractions lie in the range (9 — 
60) X 10-6 1]. 

The all-neutral-meson final states studied here are de- 
scribed theoretically by suppressed amplitudes, with pre- 
dicted branchingfractions less than a few per million by 
most estimates UlllSllllIllil^- By bringing the 
experimental sensitivity down to this level we can test 
and constrain the models. In particular, these branching 
fractions or limits bear on the accuracy with which CP- 
violating asymmetry measurements can be interpreted. 

Theoretical approaches include those based on flavor 
SU(3) relations among many modes |^ effective 
Hamiltonians with factorization and specific i?-to-light- 
meson form factors 0, perturbative QCD i, and QCD 
factorization 0- The decays to combinations of 77^'^ and 
Lo involve color-suppressed tree, CKM-suppressed pen- 
guin, and flavor-singlet penguin amplitudes, while only 
the last of these contributes to those with a single 4> 
meson. The B'^ — > (fxj) decay is a pure penguin anni- 
hilation process with an expected branching fraction of 
order 10"^ in the Standard Model this mode would 
therefore be particularly sensitive to physics beyond the 
Standard Model. 

In the time evolution of B^ v'^s ^^"^ ~* 'P^s ^ 
sinusoidal term arises from interference between decays 
with and without mixing. The coefficient S of this term 
is related to the CKM phase P = aig{-VcdV;^/VtdV*i,) if 
these decays are dominated by the single amplitude ex- 
pected in the Standard Model. Additional higher-order 
amplitudes with different weak phases would lead to de- 
viations AS between the value measured in these rare 
modes and the precise determination in the more copi- 
ous charmonium K'^ decays. Flavor SU(3) 0, 0| relates 
the strength of such additional amplitudes to the decay 
rates of two-body i?° decays to final states containing 7r°, 
r], and rj' . The 77*-'^ combinations reported here provide 
the strongest constraints. 

The results presented here are based on data collected 
with the BABAR detector at the PEP-II asymmetric 
e+e" collider (TJ] located at the Stanford Linear Acceler- 
ator Center. An integrated luminosity of 81.9 fb""'^, cor- 



responding to = 88.9 ± 1.0 million BB pairs, was 

recorded at the T{AS) resonance (center-of-mass energy 
y/s = 10.58 GeV). A 9.6 fb~^ off-resonance data sample, 
with a center-of-mass energy 40 MeV below the T{AS) 
resonance, is used to study background contributions re- 
sulting from e~^e~ —^qq{q = u, d, s, or c) contimmm 
events. 

Charged particles from e+e^ interactions are detected, 
and their momenta measured, by a combination of a ver- 
tex tracker consisting of five layers of double-sided silicon 
microstrip detectors, and a 40-layer central drift cham- 
ber, both operating in the 1.5-T magnetic field of a su- 
perconducting solenoid. We identify photons and elec- 
trons using a CsI(Tl) electromagnetic calorimeter. Fur- 
ther charged-particle identification is provided by the av- 
erage energy loss (dE/dx) in the tracking devices and by 
an internally reflecting ring-imaging Cherenkov detector 
(DIRC) covering the central region. 

The event selection criteria have been established with 
studies of off-resonance data and simulated Monte Carlo 
(MC) events of the target decay modes, BB, and con- 
tinuum. We select 77, 77', uj, and (j) candidates through the 
decays 77 77 (77^-y), 77 ^ tt+tt^-k^ (^tt), v' ^ 777r+7r" 

with 77 ^ 77 (77;,^^), 77' p^-i {ij'p^), UJ ^ TT+TT-TT°, 

and (j) K+K". The photon energy must be 
greater than 50 MeV for tt'^ and 77 candidates, and 
greater than 200 MeV in 77' pj. We make the fol- 
lowing requirements on the invariant mass (in MeV): 
490 < 771^.., < 600 for 

^77 J 120 < rriyy < 150 for tt , 
510 < m-^T: < 1070 for p°, 520 < ttitttttt < 570 for ?737r, 
910 < {m,f^Tr,mpy) < 1000 for rj' , 735 < 77iT^7r < 825 
for ui, and 1009 < ttix+k- < 1029 for (j>. We make re- 
quirements on DIRC measurements and dE/dx to iden- 
tify pions and kaons. Secondary tracks in 773^, 77', and uj 
candidates must be identified as pions, and in (j) candi- 
dates as kaons. 

A i?-meson candidate is characterized kinemat- 
ically by the energy-substituted mass 77zes = 
[(|s + Pq • pB)^/i?Q — p^] 2 and energy difference 
AE = Eg — ^^/s, where the subscripts and B refer to 
the initial T{AS) and to the B candidate, respectively, 
and the asterisk denotes the T{AS) rest frame. 

Backgrounds arise primarily from random track com- 
binations in e+e^ qq events. We reject these by using 
the angle 6t between the thrust axis of the B candidate 
in the T(4S') frame and that of the rest of the event. 
The distribution of |cos0t| is sharply peaked near 1.0 
for combinations drawn from jet-like qq pairs, and is 
nearly uniform for T{4S) BB events. We require 



5 



I cos^tI < 0.9. To discriminate against r-pair and two- 
photon backgrounds we require the event to contain at 
least the number of charged tracks in the decay mode 
plus one. For rj-^^rj^^ we require at least 3 charged tracks 
in the event. 

The decay mode (p<j) is very clean. Resolutions 

on Toes and AE are 3.0 MeV and 13.1 MeV, respectively. 
We define the signal region with cuts of ±3(t in AE and 
±4cr in Toes- The number of B^ — *■ <p(f) candidates in this 
signal region is 4.0j^f g. The only source of background 
is the continuum, estimated with on-resonance data side- 
bands to contribute 2.7 ± 0.4 events. 

We obtain yields in all other decay modes from un- 
binned extended maximum-likelihood (ML) fits. The 
principal input observables are AE and toes- Where 
relevant, the invariant masses TOi-os of the intermediate 
resonances, a Fisher discriminant and angular vari- 
ables H are used. For ij^^ , H,j is defined as the cosine of 
the angle between the direction of a daughter 7 and the 
flight direction of the rj relative to its parent in the 77 rest 
frame; for 7]'^^, Hp is the cosine of the angle between the 
direction of a p daughter and the flight direction of the rj' 
in the p rest frame; for lu, Ti.^ is the cosine of the angle in 
the LU rest frame between the normal to the lo decay plane 
and the i?° flight direction. The Fisher discriminant T 
combines four variables: the angles with respect to the 
beam axis of the B momentum and B thrust axis (in the 
T{4:S) frame), and the zeroth and second angular mo- 
ments Lo,2 of the energy flow about the B'^ thrust axis. 
The moments are defined by Lj —J2iPi^ l<^os 6i\'' , where 
9i is the angle with respect to the B thrust axis of track 
or neutral cluster i, pi is its momentum, and the sum 
excludes the B candidate. Further cuts on discriminat- 
ing variables and the set of probability density functions 
(PDF) used in ML fits, specific to each decay mode, are 
determined on the basis of studies with MC samples. For 
V-yiV'p'y requirement |7i^| < 0.86 is used to reduce sig- 
nificantly the background from the decay B'^ — > K*j. In 
other decays containing rj-y^ we require |7i,, | < 0.9 to re- 
move random combinations with soft photons. In rj^-yu; 
we apply a cut on the maximum 7 energy in the center 
of mass system (< 2.4 GeV) to suppress cross-feed from 
other BB decays with energetic photons, and a 7r° veto 
to suppress potential cross-feed from W7r°. 

We estimate BB backgrounds using simulated samples 
of B decays. The branching fractions in the simulation 
are based on measured values or theoretical predictions. 
The estimated BB background is negligible. 

For each event i and hypothesis j (signal or continuum 
background), the likelihood function is 



C 



Nl 



N 

n 



(1) 



where N is the number of input events, rij is the num- 
ber of events for hypothesis j and Vj{'x.i) the corre- 



sponding PDF, evaluated with the observables of the 
ith event. Since the correlations among the observables 
in the data are small, we take each V as the prod- 
uct of the PDFs for the separate variables. We deter- 
mine the PDF parameters from simulation for the sig- 
nal and from sideband data (5.20 < toes < 5.27 GeV; 
0.1 < |A£;| < 0.2 GeV) for continuum background. We 
float some of the continuum PDF parameters in the max- 
imum likelihood fit. We parameterize each of the func- 
tions 7^sig(wEs), 'PsigiAE), Vj{T), and the peaking 
components of Vj{m^es) with either a Gaussian, the sum 
of two Gaussian distributions, or an asymmetric Gaus- 
sian function as required to describe the distribution. 
Slowly varying distributions (mass, energy for combina- 
toric background and angular variables) are represented 
by linear or quadratic dependencies. The combinatoric 
background in toes is described by the ARGUS function 
— x"^ exp [—^(1 — x^)] , with x = 2toes/\/s and pa- 
rameter ^. Large control samples of B decays to charmed 
final states of similar topology are used to verify the sim- 
ulated resolutions in AE and toes- Where the control 
data samples reveal differences from MC in mass or en- 
ergy resolution, we shift or scale the resolution used in 
the likelihood fits. The bias in the fit is determined from 
a large set of simulated experiments, each one with the 
same number of qq and signal events as in data. 

In Table m we show the measured yield, the efficiency, 
and the product of daughter branching fractions for each 
decay mode. The efficiency is calculated as the ratio of 
the numbers of signal MC events entering into the ML fit 
to the total generated. We compute the branching frac- 
tions from the fitted signal event yields, reconstruction 
efficiency, daughter branching fractions, and the number 
of produced B mesons, assuming equal production rates 
of charged and neutral B pairs. We correct the yield 
for any bias measured with the simulations. We com- 
bine results from different channels by adding the values 
of — 21n£, taking account of the correlated and uncor- 
related systematic errors. We report the statistical sig- 
nificance and the branching fractions for the individual 
decay channels, and for the combined measurements also 
the 90% C.L. upper limits. 

The statistical error on the signal yield is taken as the 
change in the central value when the quantity — 21n>C 
increases by one unit from its minimum value. The sig- 
nificance is taken as the square root of the difference be- 
tween the value of —2 In/! (with systematic uncertainties 
included) for zero signal and the value at its minimum. 
The 90% confidence level (C.L.) upper limit is taken to 
be the branching fraction below which lies 90% of the 
total of the likelihood integral in the positive branching 
fraction region. For the B^ (jxp decay mode the 90% 
C.L. upper limit is calculated with the Feldman- Cousins 
method [l3 |. 

In Fig. n we show projections onto toes and AE 
in the analysis of the decays t]Ijj. The his- 
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TABLE L Signal yield (before fit bias correction), detection efficiency e, daughter branching fraction product, significance 
(including systematic errors), measured branching fraction B, and 90% C.L. upper limits (UL) from this and previous work. 



Mode 


Yield 


6(%) 


ne. (%) 


S(a) 


B(10-®) 


This UL (10"^) 


Previous UL (lO"**) [Ij 




'•^-5.9 


21.6 


15.5 


0.0 


^ .4-9 '^ 






^77 '?37r 
?737r 77371 


^ ^+6.8 

— n 1 


16.9 
12.3 


17.9 
5.1 


0.1 
0.0 


0.4t^-J 
-0.4«-J 






■nv 








0.0 


-o.g^i:^ ± 0.7 


< 2.8 


< 18 


'777'7r77r7r 
VllVp'Y 


-7 1 +3.V 

0.6^^3 


21.5 
20.2 


6.9 
11.6 


0.0 
0.2 


/i+2-y 
-2.4Ii g 






1 

V^TT'lrj-K-K 


1 0+7-7 


13.7 

1 s 
10.0 


4.0 


1.0 
u.o 


8.0 7 3 

9 c+S-l 






vv' 








0.3 


0.6^17 ± 1.1 


< 4.6 


< 27 


1 1 

W W 

'lr}7V7r /p7 

rj'v' 


n + 2. b 

0-3_i.5 

*-U_6.2 


14.1 
1 9 7 


3.1 
10.2 


0.1 
0.4 


0-2_4.o 

9+6.4 

1.7+^;? ± 0.6 


< 10 


< 47 




24.2^7.1 

r, 9+9.4 
^•^-8.2 


18.1 

1 9 Q 


35.1 
9n 1 


5.1 
u.o 


A A-\~i.'o 

4-4Ii.3 
9+''-i 






•quj 








4.3 


4.0+^ 2 ± 0.4 


< 6.2 


< 12 




Q+4.9 
'^•''-3.4 

^■^-i.O 


14.0 


10. D 

9fi 


n n 
U.U 

n 9 


1 0+2. b 
4+^ '-* 






•q' uj 








0.0 


-o.2tl:i ± 0.4 


< 2.8 


< 60 


'777 


-f 0.1^^9 

-2.ot?:« 


29.7 
20.9 


19.4 
11.1 


0.0 
0.0 


-2 0+^-" 
^•'-'-0.7 

-0 9+^ * 














0.0 


-l.4^°:^±o.2 


< 1.0 


< 9 


»7p70 




23.2 
22.0 


8.6 
14.5 


0.1 
1.2 


0.3+1-7 

2.8t^;« 














0.8 


l.sti.s ± 0.4 


< 4.5 


< 31 


0</. 


1 q+3.2 

J^"J-1.9 


19.9 


24.2 


0.3 


0.3+^;^ ± 0.1 


< 1.5 


< 12 



their errors and by embedding MC signal events inside 
background events simulated from PDFs. The uncer- 
tainty on Ng^ is 1.1%. PubHshed data provide the 
uncertainties in the i3-daughter branching fractions (1- 
4%). Other sources of systematic errors are track recon- 
struction efficiency (1-3%) and neutral reconstruction ef- 
ficiency (5-10%). The validity of the fit procedure and 
PDF parameterization, including the effects of unmod- 
eled correlations among observables, is checked with sim- 
ulated experiments. The value of the likelihood function 
found in data is consistent with the likelihood distribu- 
tion found in simulated experiments. 

In the cfxj) decay mode the total systematic error 

is 7.6%, which we obtain by adding in quadrature the 
errors due to the different selection cuts, branching frac- 
tions of daughters, i?° production, and statistics of the 
Monte Carlo samples. 

In Grossman et al. 9], AS* = S-sm2f3 for 5° ^ tj'K^ 
is proportional (Eq. 10) to the absolute value of a param- 
eter £,n'Ks defined in their Eq. 8. A bound |Ct7'ifsl < 0-36 
is extracted via Eq. 18 from previously measured 
branching ratios to two-body combinations of tt*', ?/, and 
rj' . The present data improve this limit: |^,,'Ks| < 0.17. 

In conclusion, we have searched for eight B^ decays 
to charmless isoscalar meson pairs. We obtain evidence 
for B^ — > rycj, with a branching fraction B{B^ —> rju) = 




AE (GeV) 



FIG. 1: Projections of the candidate mES and AE for 
B° -qui. Points with errors represent data, shaded his- 
tograms the ris^uj subset, solid curves the full fit 
functions, and dashed curves the background functions. These 
plots are made with cuts on probability ratio and thus do not 
show all events in data samples. 



tograms show the data after a cut on the probability ratio 
'^sig/ ("Psig+T^bkg), where Vsig and Vhkg are the signal and 
the continuum background PDFs. The curve represents 
a projection of the PDF obtained from a fit in which the 
plotted variable was removed. 

The main sources of systematic errors include uncer- 
tainty in PDF parameterization (1-2 events) and ML fit 
bias (0.5-2 events). We estimate these errors with sim- 
ulated experiments by varying PDF parameters within 
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(4.0t};2 ± 0-4) X 10"^ with 4.3 a significance. For the 
other modes our results represent substantial improve- 
ments on the previous upper limits 0. 
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